Challenges in modelling isoprene and monoterpene emission dynamics of Arctic plants: a case study from a subarctic tundra heath

The Arctic is warming at twice the global average speed, and the warming-induced increases in biogenic volatile organic compounds (BVOCs) emissions from Arctic plants are expected to be drastic. The current global models’ estimations of minimal BVOC emissions from the Arctic are based on very few observations and have been challenged increasingly by field data. This study applied a dynamic ecosystem model, LPJ-GUESS, as a platform to investigate short-term and long-term BVOC emission responses to Arctic climate warming. Field observations in a subarctic tundra heath with long-term (13-year) warming treatments were extensively used for parameterizing and evaluating BVOC-related processes (photosynthesis, emission responses to temperature and vegetation composition). We propose an adjusted temperature (T ) response curve for Arctic plants with much stronger T sensitivity than the commonly used algorithms for large-scale modelling. The simulated emission responses to 2 C warming between the adjusted and original T response curves were evaluated against the observed warming responses (WRs) at short-term scales. Moreover, the model responses to warming by 4 and 8 C were also investigated as a sensitivity test. The model showed reasonable agreement to the observed vegetation CO2 fluxes in the main growing season as well as day-to-day variability of isoprene and monoterpene emissions. The observed relatively high WRs were better captured by the adjusted T response curve than by the common one. During 1999– 2012, the modelled annual mean isoprene and monoterpene emissions were 20 and 8 mg C m−2 yr−1, with an increase by 55 and 57 % for 2 C summertime warming, respectively. Warming by 4 and 8 C for the same period further elevated isoprene emission for all years, but the impacts on monoterpene emissions levelled off during the last few years. At hour-day scale, the WRs seem to be strongly impacted by canopy air T , while at the day–year scale, the WRs are a combined effect of plant functional type (PFT) dynamics and instantaneous BVOC responses to warming. The identified challenges in estimating Arctic BVOC emissions are (1) correct leaf T estimation, (2) PFT parameterization accounting for plant emission features as well as physiological responses to warming, and (3) representation of long-term vegetation changes in the past and the future.